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ABSTRACT 


In  this  paper,  we  study  the  asymptotic  distributions  of  the  functions 
of  the  occurrence/exposure  rates  of  several  groups  of  patients  as  well  as 
Berry-Esseen  bound  on  the  distribution  function  of  the  occurrence/exposure 
rate.  Asymptotic  distributions  of  functions  of  the  simple  risk  rates  are 
also  derived.  The  results  are  useful  in  not  only  medical  research  but  also 
in  the  area  of  reliability. 

Key  words  and  phrases:  Asymptotic  distributions,  medical  research, 

occurrence/exposure  rate,  reliability,  risk  rate. 


1 .  INTRODUCTION 


In  medical  studies,  it  is  of  interest  to  study  the  association  between 
the  occurrence  of  certain  diseases  and  the  exposure  factors.  Various  measures 
of  risk  of  a  disease  are  considered  (e.g.,  Breslow  and  Day  (1980),  Howe  (1985)) 
in  the  literature.  One  such  measure  is  the  ratio  of  the  number  of  patients 
died  to  the  total  number  of  individuals  observed  in  a  fixed  time  period.  Using 
this  measure,  various  authors  have  studied  some  of  the  statistical  problems 
connected  with  the  risk  rate.  Another  measure  used  in  the  literature  for  the 
risk  is  the  ratio  of  the  number  of  persons  died  to  the  total  number  of  years 
exposed  to  risk.  For  surveys  of  some  developments  on  the  theory  of  occurrence/ 
exposure  rates,  the  reader  is  referred  to  Berry  (198?)  and  Hoem  (1976).  The 
main  object  of  this  paper  is  to  study  some  problems  connected  with  the  occurrence/ 
exposure  measure.  Some  results  are  also  obtained  on  risk  rates. 

Suppose  an  experiment  is  conducted  for  a  fixed  period  of  time  T  and  n 
patients  are  observed  during  this  period.  Also,  let  denote  the  total  time 
i-th  patient  is  exposed  to  risk.  Then,  the  risk  measure  considered  in  this 
paper  is 


Rn  =  VUn 


where  Un  =  Yj  ♦...+  \  =  lx  Z„. 


if  Xi  <  T 


if  Xi  >  T 

if  Xi  <  T 
if  Xi  >  T. 
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The  denominator  in  (1.1)  is  known  as  person-years. 


In  Section  2  of  this  paper,  we  establish  asymptotic  normality  of  a 

function  of  R  .  In  Section  3,  we  establish  the  Berry-Esseen  bound 
n 

on  the  distribution  of  Rn<  This  bound  is  quite  useful  since  it  gives  an 
upper  bound  on  the  absolute  value  of  the  difference  between  the  distribu¬ 
tion  functions  of  R  and  the  normal  variable  with  mean  zero  and  variance 
n 

one.  The  bound  is  of  order  c//n  where  c  is  a  constant  and  n  is  the  sample 

size.  The  asymptotic  distributions  of  the  ratios  of  the  measures  in  several 

groups  are  given  in  Section  4.  In  Section  5,  we  consider  the  measure  Vn/n 

and  give  results  analogous  to  those  given  in  Sections  3-4  for  the  measure 

Rn .  The  results  of  this  paper  are  useful  not  only  in  medical  research  but 

also  in  the  area  of  reliability.  For  example,  consider  the  situation  when 

n  items  of  an  equipment  are  under  test  for  performance  under  stress  over  a 

period  of  time  T.  A  measure  of  reliability  of  the  equipment  is  the  ratio 

of  the  number  of  items  which  did  not  fail  to  the  total  number  of  items  under 

test  during  the  period  of  time  T.  It  is  also  of  interest  to  find  the  ratio 

of  the  number  of  items  which  did  not  fail  to  X.  +  ...  +  X  where  X.  denotes 

1  n  i 

the  duration  of  the  time  i-th  item  is  under  test. 


\  wT  Of  lAhi 


2.  ASYMPTOTIC  NORMALITY  OF  THE  OCCURRENCE/EXPOSURE  RATE 

Let  p  =  P[X . >T]  =  1-q.  If  p  =  1,  then  Rn  =  0  whereas  Rn  =  n/ 

(Xj  +...+  Xn)  when  p  =  0.  Both  of  the  above  cases  are  simple  and  so  we 
deal  the  case  when  pe(0,l). 

Using  strong  law  of  large  numbers  for  i.i.d.  sequence  ,  we  have 
(Vn/n)  -*■  q  almost  surely  (a.s.)  and 

n  Un  -  M  =  EYX  *  E(X1)ICX1<_T]  +  Tp  a.s. 

as  n  -*■  ».  Hence,  Rn  -*■  q/u  a.s.  Now,  let  *  uZ^  -  qY^  ,  r=q/u 
and 

5n  =  ^  (Rn'r) 

“HD"  .1  wi3-  (2.1 

n  / n  i=l 

Here  {W. }  is  a  sequence  of  bounded  i.i.d.  random  variables  with  mean  zero. 

n 

So,  by  central  limit  theorem,  we  observe  that  £  W.//rT  is  asymptotically 

i=l  1 

2 

distributed  as  normal  with  mean  zero  and  variance  a  ,  where 

o2  =  E(W 2)  =  EjuZj-qYj)2.  (2.2 

1  2 

Since  -  uUp  ->  u  a.s.,  we  obtain  that  Cn  is  asymptotically  distributed  as 

2  4 

normal  with  mean  zero  and  variance  a  /u  . 
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Now,  let  f(-)  denote  a  function  which  is  continuously  differentiable 
for  two  times  around  r,  say  in  (r-6  ,  r+6),  6  >  0.  By  Taylor's  expansion 


if  |Rn-r|  <  6/2,  we  obtain 


(f(Rnl-f(r)) 


where  cn  is  a  number  between  r  and  Rn*  Because  f"  is  bounded  in  the 
interval  (r-6/2  ,  r+6/2),  tends  to  a  normal  variable  in  distribution 

and  P(  I  Rn-i“  |>6/2)  -*  0,  and  we  have  the  following  theorem. 

THEOREM  1 


Under  the  condition  mentioned  above, 

/T  (f(Rn)-f(r)^  N  (o,  (f'(r|2a2/u4^. 

In  practice,  the  asymptotic  variance  of  /n”  ^f(Rn)-f(ry  is  unknown. 
In  such  situations,  we  use  the  following  approximate  confidence  interval  on 


^  (f(Rn)-f(r))|  l  doa(f) 


where  a(f)  can  be  taken  as 


which  is  a  consistent  estimate  of  |f'(r)|a/u^  and  dQ  is  the  upper  100a% 
point  of  the  normal  distribution  with  mean  zero  and  variance  one. 


3.  BERRY-ESSEEN  BOUND  FOR  THE  DISTRIBUTION  OF  THE  OCCURRENCE/EXPOSURE  RATE 


Let 


i  n 

l 

/TT  o  i  =1 


W... 


Then,  according  to  the  result  proved  in  previous  section,  np  is 

asymptotically  distributed  as  normal  with  mean  zero  and  variance  one. 

Let  F  denote  the  distribution  function  of  n  and  4>  that  of  the 
n  n 

standard  normal.  In  this  section,  we  shall  prove  the  following. 


THEOREM  2 


There  exists  a  constant  c  such  that 

UF  -  <5>U  =  sudIF(x)  -  $(x)l  <c/SrT  .  (3.1) 

n  x  n 

where  k  is  the  standard  normal  distribution  function.  In  the  sequel, 
we  need  the  following  lemma. 

LEMMA  1 


Let  { X n ,  Yn,  Zn>  be  a  sequence  of  random  vectors  with  relation 
s  Yn  +  Zn  and  let  Fn,  Gn  denote  the  distribution  functions  of  Xn 

and  Yn  respectively.  If  there  exist  constants  ci ,  i  =  1,2,3,  such  that 

1)  I  !Gn  -  *1 1  1  Cj/Zn". 

2)  P(|Zn|  >  c2//TT)  <  c3//F  , 
then  there  exists  a  constant  c4  such  that 


>va; 

m 


v.‘ 

m 

V.-- 

■ ■ 


P7K 

.yjw 

v\^. 


1 1  Fn  "  $il  1  c4//rT . 


For  a  proof  of  the  above  lemma,  the  reader  is  referred  to  X.R.  Chen  (1981). 


Now,  we  turn  to  prove  Theorem  2.  Let 


Tn  -  1  +  W  ^ 


n_  =  S  +  A*  +  a"  +  a!"  +  a  ' l* '  , 
n  n  n  n  n  n  * 


(3.2) 


(3.3) 


where 


i  n 

h-  i 

no  i=l  1 


sn  ■  — 

n  /rTc 


An  3/2  I  W.(l-Y./u), 
n  nJ/^o  l<i?j<m  J 


An  =  ~TTt  l  W.(l-Y./u), 
n  o  i  -l 

An'  =  n3^  ^1  Wj^"Yi/u^ 


A  II  I  I.  A  C 

A_  - 

n  n  n* 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


*nd  m  =  n  -  /n“,  the  summation  runs  over  all  possible  values  of  i  and  j  such 
that  l^i<.n,  m+l_<  j  <n ,  i£j  or  l_<j£n ,  m+l<i<n,  i  j4  j . 

At  first,  we  see  that 

P(l^  I  >  T=)  1  nE(An)2 

/n 

*  o_2n"2  I1[EW?(l-Yi/u)2  +  2EW.Wj(l-Y./u)(l-Yj/u)] 

<  3o'2n~1/2  EW2(1-Y2/u)2 

<  c/.'T,  (3.9) 


where  and  in  the  sequel  c  denotes  positive  constant  but  may  take  different 
value  at  each  appearance.  Also,,  for  any  c  >_  a-1  EW^l-Yj /u)  |  +  1^,  we  have 

P^|A„I  >  c//n 
=  P^l  ^Wi(l-Yi/u)|  >  Can^ 

<  p(l  f  w. (l-Y./u)  -  EMjd-Yj/uJI  >  n) 

<  h_1Var  ^Wjd-Yyu)^  <  c//F  .  (3.10) 

We  now  estimate  A1'".  Define  the  event 

E"=  (|  ",ii(1'Vu>l 

L>j  Hoeffding  inequality  (see  Hoeffding  (1963)),  we  have 

P(En)  <  2  exp{-  2n(l/2T)2}.  (3.11) 

Let  E^  denote  the  complement  of  the  event  En>  When  E^  is  true,  we 
have 

Thus 

P{|A'’M  >  1//TT) 

=  Ansnl  >  1//F) 
i  P<E„>  +  pCv  IVnl  i 
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te 
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-  P<En'  *  P(2(n  1-I1(1"V1/u0  Z[  TFo  1  “l1-1^) 

i  PtEn)  *  P  (|i  I  (1-V i/u)|  >-L  »"*Y  p/|-A.  JW,|  »„M(3.12) 
\  i=l  /2  /  \/rT a  i=l  / 


By  Hoeffding  inequality,  we  get 

li 


and 


('*  j1(1-vi/u)l-^"'3/8) 


1  2  exp{-  2  n(— — n  -3/8)2} 

/Tt 


l0/^  (3.13) 

1  2  exp{-  n"1/4))2>  <  c//rT  .  (3.14) 

From  (3.11)  -  (3.14),  it  follows  that 

P^|a“"|  1  l//rT^  <  c//7T 

Applying  lemma  3.1,to  prove  Theorem  2,  we  only  need  to  prove  that 

I ! Gn  -  $|[  <  c//7T  ,  (3.15) 

where  G  denotes  the  distribution  function  of  T  =  S  +  A1  ,  and  S  , 
n  n  n  n  ’  n 

a^  were  defined  in  (3.4),  (3.5). 
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V.v\v..-,s 


■.  s 


(The  proof  of  (3.18)  can  be  found  in  Chapter  5  of  Petrov's  book  (1975)). 
Therefore  there  exists  a  constant  5^  >  0,  such  that  for  any  |t|  £  <Sj , 

I v (t ) |  <  exp{-  ^  t2}.  (3.19) 

Hence  for  jt|  £  <S^/n”,  we  have 

|Ee1tSn2|  <  |v(t//n-)l[/ir] 

£  exp{-  t2/4/n~}  (3.20) 

and 

|E  exp{it(Sn-(a1+^)//F  )}| 

=  |v(t//rT) |n_2  <  exp{-  t2(n-2)/4n} 

<  exp{-  t2/5>  for  large  n.  (3.21) 

Note 

E(An)2  =  "T7  m(m“l)[Eb?2  +  2Eb12b21] 
no 

<  c/n. 

Hence  from  (3.17),  we  get  for  | t |  £  6^/17  , 

2  itS 

|  E(A^)29ne  n|  £  c(t2/n)  exp{-  t2/4/F  }.  (3.22) 

Now  write 

gn(t)  =  Eb12  exp{it(a1+a2  )//rT  )  . 


Since  Eb^  =  0.  =  ^*3l2a2  =  we  *iave 

l9n(t)|  -  In  Edbi2l  )(ai+a2^ 

<  ct2/n. 


(3.23) 


By  (3.21),  (3.23),  we  have  for  |t|  <.  5j  /n~  , 

'EAne  n|  1  l9n(t)|  |E  exp{it(Sn-(a1fa2)//TT)}| 


n  a 


<_  Ct2n-1^2  exp{-  t2/5> .  (3.2 

From  (3.16),  (3.22)  and  (3.24),  we  get 

|f„(t)-fn(t)|  <  c^-V*^5  e-t2/4/T).  (3.2 

/n 

By  lemma  1  in  Chapter  5  of  Petrov  (1975),  we  have  for  |t|  <  62/lT,  >  0 

|fn(t)-e-t2/2|  <cit^e-t2/8.  (3.2 

n  /F 

Thus  (3.25)  and  (3.26)  yield  for  |t|  <  5  /IT 

|fn(t)-e”t^/2|  <  c[(i^  e’t2/8  +  ^  e_t2/4/7r]  (3*2 

/~n 

where  s  =  min(61,60)  >  0.  From  (3.27),  it  follows  that 
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Here  the  estimate  of  the  last  integral  can  be  obtained  by  making  variable 
transformation  u  =  tn-1^4.  Then  using  Berry-Esseen's  basic  inequality, 
we  prove  (3.15).  This  completes  the  proof  of  Theorem  2. 


14 


4.  ASYMPTOTIC  JOINT  DISTRIBUTION  OF  FUNCTIONS  OF  OCCURRENCE/EXPOSURE  RATES 
Let  x{J) . xW  ,  (j*=l,2 . s)  be  a  sample  drawn  from  the  j-th 

J 

population  where  x|J^  denotes  the  observation  on  i-th  individual  in  j'-th 
population.  Also,  let 


rU)  .  1 


if  <  T 


otherwise 


,(j)  _ 


if  <  T 


otherwise. 


for  j  =  l,2,...,s  and  i»l,2,...  .n^.  Now,  let 


R(j)  s  v( j )/L|( j ) 

"j  nj  "j 


for  j=l,2,.. .  ,s,  where 


uU)  =  nf  Y(j)  v(j)  =  Z(J). 

nj  i=l  1  nj  i=l  1 


(4.1) 


(4.2) 


We  know  that 


Rn.  rj  a,s*  j  =  l»2,...,s  . 

J  J 


(4.3) 


Let  f(xj,X2»...,xs)  be  a  function  which  is  continuously  differentiable 
for  two  times  in  a  neighborhood  of  (rj,...,r$).  Suppose  that 


n/n .  ■+  a.  <®  ,  as  n  , 

J  J 


(4.4) 


where  n  =  nj  +,...,+  n  j . 


,vy. 

x-x 

NCyV. 


■  w  V  -  •  -  —  -  - 


/7T  (f(R^»  Rn^*“**Rn^  “  f^rl’  r2****’rs^) 

-  y  a  JjL  E(J)  +  y  y  J  IlL  a  r(jMk) 

A  V"jS  j-1  k-1  njnk  jkS  V 


(A. 5) 


when  R^,  R^ •  »R„S^  falls  in  the  neighborhood  of  ( . r$)  in 


which  f  is  differentiable.  Here 


3f(x1,...,xs)| 


ajk  =  3xTxT 


(Xj»»..»Xj)  -  (rj,...  ,fj)  $  j  -  1.2 . s. 


(x1’-*’»xs)  s  . V  .  j,k  =  1.2. . s 


and  (tj,...,ts)  is  some  point  on  the  linear  section  joining  R^j}. . .  .R^ 

and  (rj,...,rs).  Let  B  be  a  non-trivial  closed  ball  with  center 
(r^,...,rs)  which  is  contained  in  that  neighborhood  of  (r^ . r$).  Then, 


P((Rij)--”Rn*))  *  B)  -  °* 


Since 


|a.kl  <M 


for  all  j  ,k  =  1,2,  ...,s  and  some  M  when  (R^ ,. . .  ,R^sb  e  B,  we  obtain 

nl  ns 

pf  I  y  ya  JUI  E(J)g4k)|  „  \ 

A-i  a  JkS \  S  k  -1  - ; 

<  »)  *  jj  j,  p(i4j>^)i  °- 


r.  r.v. r.  •</ 
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Hence 


/F 


(f(R^1} . R^})  -  f(r j . r$j)  -  N(0  ,  o2)  (4.6) 


where 


’?  ■  1  aj  *j  °Wi  ■ 


u,  =  E  Y]  , 


,(j) 

rl 


q.  =  P(x{j)  <  T) 


and 


Oj  =  E(wJj))2  =  EfUjZ^  -  qjY{j))2. 

.  Also,  using  the  same  approach  to  prove  Theorem  2,  we  can  establish 
the  Berry-Esseen  Bound  for  the  distribution  of  ^f(R^ ,R^)  - 

,rs)J  .  The  details  are  omitted  here. 

An  important  special  case  for  f  is  f(Xj,x2)  *  x^/x^.  In  this  case, 

fCR^.R^)  =  R^Vr^  is  called  the  ratio  of  occurrence/exposure  rates. 


nl  "2 


nl  "2 


rO)/r(2)  1<s  (ignored  by  ,  and  we  have 


nl  n2 


where 


and 


(RRn  -  ryr2)  -  N(0,  aZ) 


2  -  fl2,  2 ,..4  .  .2,  Z..A 
cj  —  ^  a2^2a2' U2 


=  1/^  »  a2  =  “  ^i/r2 


(4.7) 


^--;X 

•  •.*  .  ^ . 
V.  W 

I*"* 

»  .  >r.  *• 

<VvS 


Note  that  R^'  may  be  zero.  ^However,  P(R^  =  0)  =  p2  -*•  o,  as  n.-*- 

Any  way,  the  definition  of  RRn  for  rJ,2)  =  0  does  not  affect  the  limiting 

result  for  the  distribution  of  RRr.  However,  for  small  sample  problem,  we 
have  to  make  an  explicit  distribution  of  RR  =  ®  when  R^2'  =  0.  Define 

RRn  =  *  w^en  Rnj^  =  Rn2^  =  ®  anc*  RRn  =  “*  ^ow  commen  density  of 


xjJ^,...,X^  be  given  by 
J 


9j(x)  " 


0j.  exp{-  OjX)  if  x  >  0 


otherwise 


and  let  Pj  =  P[xjj)  >  T]  for  j  =  1,2.  We  have 

A 


✓V.  II,  n7 

p(RRn=o)  =  p^Ci-Pa2) 

A  n  n 

P(RRn=l)  =  p/pg2 

^  n,  n, 

P(RRn=«)  =  P2 d-Pi  ) 


J 


(4.8) 


It  is  known  (see  Beyer,  Keiding  and  Simonsen  (1976))  that  has  an  atom  at 

n .  J 

the  orgin  with  a  mass  P^J  and  a  density 

fj(X)  k-1  ^  PJJ  kotjTk’lkx"2exP{-  «j(k-(n-k)Tx)+/x) 


uk[(k-(nj-k)Tx)+/x]  I 


[0,  k/(nj-k)T] 


(x). 


(4.9) 


for  j  *  1,2,  where  I^a  ^(x)  is  1  or  0  according  as  x  is  in  [a,b]  or  not.  Henc 
the  distribution  of  RRn,  besides  the  three  atoms  given  in  (4.8),  has  a  density 
which  can  be  computed  from  the  following 

tm 


f  (x)  = 


fl(xy)f2(y)dy 


(4.10) 
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Now,  let 

Li  *  JT  (fi(Rj1j),...,R^))  -  . r$)} 

for  i  =  1,2, f,  (R^ . R^)  is  a  continuous  twice-differentiable 

1  nl  ns 

function  of  R^ ,. . .  .R^  arount*  rj,...,r$.  We  have  proved  earlier  the 

asymptotic  normality  of  Lj.  Following  the  same  lines,  it  is  easily  seen 
that  the  asymptotic  joint  distribution  of  Lj,...,Lk  is  multivariate  normal. 

But  the  asymptotic  covariance  matrix  of  Lj . Lk  is  isually  known.  We 

will  now  construct  approximate  confidence  intervals  on  f- (rj,... ,r$)  when 
the  covariance  matrix  C  =  (c^)  of  Lj,...,Lk  is  non-singular,  where 

2,4 


and 


:it  ■  i.  ai-jat-jxj°j/uj 


i*J 


3fi(xi .... ,XS)‘ 


(Xj ,. . . ,x^ )  =  (r^ ,. . .  ,r$) 


a  *  1 

In  these  situations,  let  C  be  a  consistent  estimate  of  C.  Then  L * C  L 

is  approximately  distributed  as  chi-square  with  s  degrees  of  freedom  for 
large  samples  where  L‘  *  (Lj,...,Lk).  Using  this,  we  obtain  the  following 

approximate  confidence  intervals  on  linear  combinations  of  f^(rj . r$)  , 

(i  =  1,2,... ,k): 


[/Fa,(f(R!j), 

for  all  nonnull  vectors  a:kxl  where 


,R<*})  -  HrY . rs))\  <  (9aa-Ca)1/2 


A 


fCrj,...,^)  =  (fjJrj,...  ,r$)  ,...,  fk(r1#. >rs)) 

and  g^  is  the  upper  100a%  point  of  the  chi-square  distribution  with  s 
degrees  of  freedom.  The  above  confidence  intervals  are  useful  in  constructing 


P 


.*  s  , 
.*,  r.  +' 
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simultaneous  confidence  intervals  on  various  ratios  like 


r../rs  (i  =1 ..  ,k-l)  ,  ri/rJ. ( i  <  j  =  2 . k) 


r*/r<_j_^  (i=l»2,...  ,k-l). 


r.w, .  w; 

W 

W 

a/v 

>>!• 

V’.C) 


He  can  also  construct  simultaneous  confidence  intervals  on  fj(rj»...»r  ) 
using  Bonferroni's  inequality. 


>■9 
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5.  INFERENCE  ON  SIMPLE  RISK  RATES 


In  this  section,  we  compare  the  simple  risk  rates  of  different  groups  of 

patients  who  are  observed  for  a  fixed  period  of  T  years  and  each  group  may 

be  subject  to  a  different  exposure  factor.  Here  a  simple  risk  rate  of  j-th 

population  is  defined  as  the  proportion  of  individuals  in  that  population  who 

died  during  the  period  of  observation.  In  this  section,  we  use  the  same 

notation  as  in  the  preceding  sections. 

The  sample  estimate  of  simple  risk  rate  for  j-th  population  is  V.  = 

Now,  let  f.j  (Vj ,. . .  ,Vs)  ,  (i=l  ,2 ,. . .  ,k),  be  a  continuous  twice  differentiable 
★  ★ 

function  of  Vj,...,V  around  qj . qs. 

Using  Taylor's  expansion,  we  obtain 


L1  *  JT  {f.(V*,...,V*)  -  fi  (q2  » —  ,q$)> 


s  s 


■  ji,  aij  *7=^  h  (5-'> 


where  =  /nT  (V^/n^)  -  ,  Vj  =  vj^  ,  and 

J  J 

=  8fi  .  s2fj 

3Vj  v  =  q  3Vj.8Vk 


3Vj8Vk  K  =  5 


(5.2) 


V  =  (Vj,...,V  )'  and  q  =  (q^,...,qs)'  and  a  is  some  point  on  the  linear 

★ 

section  between  q  and  V  .  As  n  -*•  ®>,  B.  is  distributed  as  normal  with 

*»  **  sJ 

mean  0  and  variance  q.p..  So,  when  n,  n,,...,n  -►  the  joint  distri- 

J  J  x  s 

★  ★ 

bution  of  Lj,...,Lk  is  multivariate  normal  with  mean  vector  0  and  covariance 
★  ★ 

matrix  C  =  (cjt)  where 


'it  =  aijatjXjqjpj  ' 


(5.3) 


vs  v-vS  \,s 


r^.  vvcvjt:?1;  v-j  *-_r; 
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Let  C  be  a  consistent  estimate  of  C  .  When  C  is  non-singular  and 


n  -*•  »,  we  can  use  the  following  approximate  simultaneous  confidence  intervals 

*  I  -  * 

for  the  linear  combinations  of  q^ . qs  by  using  the  fact  that  V  CV  is 

approximately  distributed  as  chi-square  with  s-  degrees  of  freedom 

!  a'(f(vi . -  f (qj  » —  »qs^)  I  1  (ha?'^*!^1/2  (5.4) 

where  f(q1#...  ,q$)  *  (fjfqj ,. ..  ,q$)  . fk(qj . q$)) '  and  ha  is  the 

upper  100a%  point  of  the  chi-square  distribution  with  k  degrees  of  freedom. 

ic  ★  &  ★  ★  ★ 

Some  special  cases  of  f. ,V  )  are,  1^/Vg  ,  v-j/Vi+j*  etc.  From 

the  results  given  above,  it  is  easily  seen  that  /n^+n^  ^R12  -  (qjA^)^  1S 

2 

distributed  normally  with  mean  zero  and  variance  aQ  where 

°0  =  +  X2^1c*2^2^<^2  *  ^12  =  ^1^2 
when  n1  and  n2  tend  to  infinity.  Following  similar  lines  as  in  Section  3, 
we  can  show  that 


lFV"2  -  i 


/hj+hg 


where  Fn  +n  is  the  distribution  function  of  t^+n^  o"1  Rj2  -  (qjA^)  and 

*  is  the  distribution  function  of  the  standard  normal  distribution. 

We  know  that  follows  the  binomial  distribution  B(n.  ,  q.),  j  =  1,2,.. 

*'j  J  J 

whatever  the  underlying  distributions  are.  Hence,  we  have 
/ 


(i-qi)  Hi-d-qg)  z3< 


p(r12  -  x)  -  < 


(1-U-q/1)  (l-qg)"2 

Here,  the  summation  ^ j  runs  over  all  possible  values  of  k^  and  k2  such  that 
1  i  k  _<  rij  ,  1  _<  k2  <  n2  and  ( k1/n1)  =  x(k2/n2)  and  the  term  for  k1  =  k2  =  0 
appears  only  when  x  =  1. 
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if  x  =  0, 
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